均等 分 区 问题 模型 、 算 法 及 应 用 
孔 云 峰 : 


1 河南 大 学 黄河 中 下 游 数字 地 理 技术 教育 部 重点 实验 室 ， 河 南 开封 ，475000 


摘要 : 分 区 问题 广泛 应 用 于 地 理 、 经 济 、 政 治 、 商 业 、 公 共 服 务 等 领域 。 均 等 分 区 问题 是 其 中 一 
类 问题 ， 通 常 要 求 分 区 人 口 数 量 或 任务 量 均等 、 几 何 形 状 紧 次 和 空间 连续 ， 应 用 于 选区 、 销 售 区 
和 巡逻 区 的 划分 。 本 文 针 对 均等 分 区 问题 提出 了 一 个 混合 整 型 线性 规划 模型 ， 并 设计 了 一 个 基于 
迭代 局 部 搜索 (ILS) 的 混合 算法 。 该 算法 从 三 个 方面 扩展 ILS: 群 解 搜索 、VND 搜索 及 SPP 模型 
改进 。 选 择 5 个 区 域 对 模型 和 算法 进行 测试 ， 结 果 表 明 : 数学 模型 能 够 求解 空间 单元 数量 为 324 
的 案例 ; 混合 算法 优化 性 能 优异 ， 鲁 棒 性 强 ， 计 算 效 率 较 高 。 所 提出 的 均等 分 区 问题 适用 于 政治 


分 区 等 经 典 问题 ， 在 新 冠 疫情 应 急 服务 等 领域 也 具有 应 用 潜力 。 
关键 词 : 均等 分 区 问题 ; 混合 整 型 线性 规划 ; 混合 算法 ; 应 急 服 务 


分 区 问题 是 将 特定 地 理 空间 划分 为 若 
干 区 域 , 满足 一 定 的 约束 条 件 ， 并 达到 最 优 
的 分 区 目标 , 广泛 应 用 于 地 理 、 经 济 、 政治 、 
环境 、 生 态 、 农 业 、 城 市 、 商 业 、 公 共 服 务 、 
地 图 制图 、 空 间 统 计 等 领域 。 分 区 问题 类 型 
繁多 , 涉及 分 区 数量 、 单个 或 多 个 分 区 指标 、 
静态 或 时 间 序 列 指标 、 是 否 要 求 分 区 空间 连 
续 等 。 根据 分 区 遵循 的 基本 原则 ,分 区 问题 
可 分 以 三 大 类 : () 均 质 性 分 区 , 即 分 区 内 差 
异 最 小 化 和 分 区 间 差 异 最 大 化 ， 如 自然 地 
理 、 气 候 和 农业 区 划 ; (2) 功 能 性 分 区 ， 即 分 
区 内 部 交互 作用 最 大 化 (区 内 凝聚 ) 和 分 区 间 
交互 作用 最 小 化 (区 间 分 离 )， 如 经 济 功 能 区 
和 城市 功能 区 划分 ; (3) 均 衡 性 分 区 ， 如 分 区 
人 口 数量 或 其 他 属性 均等 的 选区 、 警察 巡 远 
区 和 销售 区 划分 , 供应 与 需求 均衡 的 学 区 和 
医疗 服务 区 划分 。 自 1960s 以 来 ， 学 术 界 一 
直 在 探索 各 类 分 区 问题 求解 方法 , 相关 算法 
研究 文献 十 分 丰富 。 然 而 ， 由 于 分 区 问题 涉 


制 要 求 , 基于 人 口 统计 基本 单元 将 一 个 区 域 
划分 为 若干 选区 ， 在 选区 中 进行 政治 选举 。 
选区 划分 通常 优先 考虑 选区 人 口 均等 、 空 间 
连续 和 几何 形状 紧凑 原则 ,从 而 避免 选区 划 
分 的 不 公平 性 ， 其 他 准 测 (如 与 行政 区 吻合 
度 ， 种 族 、 社 会 、 经 济 、 政 治 属性 均衡 分 布 
等 ) 并 未 形成 共识 [1]。 销 售 区 划分 问题 是 企 
业 将 其 市 场 区 域 划 分 为 在 干 分 区 ,指派 一 个 
雇员 或 一 个 团队 负责 一 个 区 域 进 行 产 品 销 
售 或 客户 服务 ， 从 而 合理 、 高 效 地 进行 企业 
运作 管理 。 销 售 区 划分 首先 考虑 雇员 工作 量 
均等 ， 其 次 是 分 区 紧凑 和 空间 连续 ， 从 而 减 
少 旅行 开支 。 类似 地 ， 和 警察 巡逻 区 划分 要 求 
分 区 的 巡逻 任 务 量 均 等 , 并 考虑 分 区 紧凑 和 
空间 连续 ， 从 而 提高 巡逻 效率 。 

现实 生活 中 的 某 些 设施 服务 区 设计 也 
要 求 服务 区 均等 、 紧 次 和 连续 。 例 如 ， 公 安 
警 务 室 设计 。 若 每 个 警 务 室 服务 范围 人 口 均 
等 、 形 状 紧凑 和 空间 连续 ， 则 使 请 警 管辖 范 


及 领域 三 ， 问 题 类 型 多 , 分 区 原则 、 目 标 函 
数 和 约束 条 件 各 有 差异 , 求解 分 区 问题 的 模 
型 和 算法 仍 在 持续 探索 之 中 。 

选区 、 销 售 区 、 警 察 巡逻 区 划分 等 问题 
是 一 类 常见 的 均等 分 区 问题 (EDP)。 此 类 问 
题 均 要 求 分 区 空间 连续 、 分 区 人 口 数量 或 任 
务 量 均 等 、 几 何 形 状 紧凑 和 基本 空间 单元 完 
整 (不 被 分 割 )。 选 区 划分 问题 是 遵循 法 律 强 
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围 明确 、 任务 量 均 衡 且 居民 办 理 相 关 事 务 的 
旅行 距离 短 。 出现 重大 疫情 实施 封闭 小 区 措 
施 时 ， 为 保证 居民 生活 物资 配送 ， 也 可 以 将 
实施 封 轩 区 域 划 分 为 在 干 服务 区 。 均 等 、 紧 
凑 和 连续 的 服务 区 设计 , 能 够 使 得 居民 生活 
物资 保障 人 员 的 任务 量 均衡 ,上门 旅行 距离 
较 短 ， 且 服务 范围 明确 。 新 冠 病毒 全 民 核 酸 
检测 亦 需 要 分 区 服务 规划 。 例 如 ， 疾 控 部 门 
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派出 多 个 工作 组 对 某 一 区 域 进行 检测 , 每 个 
工作 组 设置 一 个 服务 站 , 计划 在 最 短 时 间 内 
完成 检测 。 每 个 工作 组 的 服务 区 要 求 任务 均 
等 、 服 务 范 围 紧凑 且 连 续 。 服 务 区 紧凑 可 使 
每 个 服务 站 服务 的 居民 空间 集中 度 高 、 到 达 
服务 站 步行 距离 最 短 , 也 能 避免 与 居住 较 远 
的 居民 产生 聚集 。 

求解 分 区 问题 的 主流 算法 包括 精确 算 
法 和 启发 式 算法 [2]。 因 分 区 问题 是 NP-Hard 
问题 [3], 数学 模型 求解 能 力 有 限 。 文献 [4] 为 
均 质 性 分 区 问题 构建 了 三 个 数学 模型 ， 然 
而 ， 针 对 49 个 单元 的 案例 ，3h 的 计算 时 间 
仍 难于 获得 最 优 解 。 文 献 [5] 构造 了 一 个 政治 
分 区 问题 数学 模型 ， 需 要 6s~53,632s 求解 
25~36 个 空间 单元 的 案例 。 为 高 效 求解 , 现 有 
模型 多 将 问题 简化 ， 如 不 考虑 空间 连续 、 不 
考虑 紧凑 目标 或 者 将 单元 分 割 [6-10]。 将 集 
合 划 分 问题 (set-partitioning problem，SPP) 模 
型 用 于 分 区 问题 ， 要 先 构造 大 量 分 区 方案 ， 
通过 模型 选择 最 优 的 分 区 组 合 [11-13]; 虽然 
集合 划分 问题 模型 相对 容易 求解 , 但 求解 结 
果 取 决 于 候选 分 区 的 数量 和 质量 。 

分 区 问题 的 元 启发 算法 很 多 , 如 基于 邻 
RE RAPAE (Greedy)、 模 拟 退 火 
(SA)、 禁 忌 搜 索 (TS) 和 动态 门槛 法 (OBA) 算 
法 [14-19]， 基于 自然 进化 的 遗传 (GA)、 演化 
(EA)、 分 散 搜索 (SS)、 智 能 蜂 群 (ABC) 算 法 
[22-28]. 然而 , 优化 性 能 较 好 的 元 启发 算法 ， 
算法 机 制 设计 较为 复杂 ， 计 算 时 间 也 偏 长 ， 
如 PEAR 算法 [28] 需 要 设置 19 个 算法 参数 ， 
在 设置 13.1 万 处 理 器 的 超级 计算 机 上 运行 
KIA 3 小 时 的 时 间 。 另 外 ， 区 域 生 长 方法 
[29][30] 和 计算 几何 方法 [31-33] 思 路 简单 ， 


对 于 一 个 地 理 区 域 ,包含 2 个 空间 单元 。 
SRE ={1,2...n} 表 示 n 个 空间 单元 ,每 一 
个 空间 单元 i 具有 属性 p:， 表 示 空 间 单 元 i 
人 口 数量 或 任务 量 ， 变 量 cj 代表 空间 单元 i 
和 j 是 否 邻 接 (有 公共 边 ); 集合 W 表示 单元 
i 的 相 邻 单元 ， 即 Ne, Vj EV, Cij = 1}; AR 
量 qj 表示 空间 单元 i 与 学 校 单 元 j 之 间 的 距 
离 。 本 文 将 均等 分 区 问题 定义 为 : 将 空间 单 
元 集合 划分 为 个 子 集 ( 即 KADE), w 
为 集合 R={7,r2...rx}， 要 求 基本 空间 单元 完 
整 、 分 区 空间 连续 、 分 区 人 口 或 任务 量 均等 、 
分 区 几何 形状 紧凑 等 。 这 一 定义 适用 于 选 
K AEK KEKEN pi. 

空间 连续 和 紧凑 的 表达 方式 是 构建 分 
区 问题 数学 模型 的 关键 。 将 地 理 空间 表达 为 
一 个 网 络 图 GV, Ey, MRE ERITI 
结 点 ， 相 邻 空 间 单元 (cy = 1) 相 连接 作为 网 
络 边 ， 可 采用 “网 络 流 ”概念 表达 分 区 空间 
连续 性 [4] [34]。 在 一 个 分 区 内 ， 每 个 单元 产 
生 一 个 单位 的 流量 ， 通 过 网 络 边 流动 ， 最 终 
汇 入 到 分 区 内 某 一 汇集 单元 。 汇集 单元 不 产 
生 流 量 , 但 汇 入 最 高 达 n-K 单位 的 流量 。 产 
流 和 汇流 单元 所 构成 的 区 域 是 一 个 空间 连 
续 的 区 域 。 SEAN; = Ulciy = Awe 
T i 相 邻 的 单元 集合 。 文献 [4] 引 入 决策 变量 
Wx 表示 单元 i 是否 分 区 上 的 汇流 单元 , 决策 
变量 玖 表示 单元 是 否 属于 分 区 k， 整 型 变 
量 fi 表示 单元 h j 间 在 分 区 中 的 流量 ， 
可 构建 线性 不 等 式 约束 保证 分 区 空间 连续 。 

分 区 几何 紧凑 度 指数 广泛 应 用 于 分 区 
几何 形状 评价 。 有 多 种 几何 紧凑 度 定 义 ， 然 
而 每 种 定义 具有 一 定 的 优势 和 局 限 [36]。 常 
见 几何 紧凑 度 指标 定义 划分 为 两 大 类 : 第 一 


但 优化 性 能 有 限 。 总 体 上 ， 现 有 分 区 问题 算 
法 往往 难以 兼顾 算法 设计 难度 、 优 化 性 能 和 
计算 效率 。 


类 从 平面 几何 角度 定义 ， 如 keR44k/TDE / 
KG，、PDPr 为 分 区 的 面积 和 周 长 ); 第 二 类 
从 空间 可 达 性 的 角度 定义 ,如 Zier, dic, (ct 为 


本 文 针 对 均等 分 区 问题 , 构造 了 一 个 混 
合 整 型 线性 规划 (MILP) 模 型 ， 满 足 分 区 均 
等 、 紧 和 和 连续 要 求 ， 设 计 了 一 个 基于 迭代 
局 部 搜索 ULS) 的 混合 算法 。 模 型 和 算法 是 、 
的 有 效 性 在 五 个 案例 区 得 到 了 验证 。 

1 问题 定义 及 数学 模型 


Min. Yiev Prev PidirXixr / Liev Pi 


区 域 的 中 心 单元 )。 本 文采 用 第 二 类 紧凑 
EEX, 即 从 每 个 单元 到 达 所 在 分 区 中 心 单 
元 的 总 距离 作为 紧凑 度 目 标 函 数 , 总 距离 越 
短 分 区 越 紧 凌 。 

将 紧凑 度 作 为 目标 , 将 分 区 均等 和 空间 
连续 作为 约束 条 件 , 构建 如 下 整 型 线性 规划 


模型 ; 


(1) 


S.T. Prev Xiz = 1,Vi E V (2) 
(1 — e)Qyx < Liev PiXik < (1 +e) Qy_, Yk EV (3) 
dev Yk = K (4) 
fix < (Nn —K) xin, Yi EV,jEN,kEV (5) 
fixn < (2 — K)xjp_,ViE Vij EN, KEV (6) 
Djen, fijr — Ljew, fjik = Xin Yİ EV,k EV, i + k (7) 
Xizx = {0,1}, vi Ee V,k EV (9) 
Yg = {0,1}, Vk EV (10) 
fin Z 0,Vi EV,j EN, KEV (11) 


模型 中 ， 目 标 函 数 (1) 是 紧凑 度 目 标 函 
数 , 分子 部 分 表示 所 有 单元 与 其 分 区 中 心 的 
加 权 距 离 , 分 母 部 分 表示 地 理 区 域内 的 总 人 


(7) 是 分 区 空间 连续 约束 条 件 ， 条 件 (5) 和 (6) 
确保 属于 同一 个 分 区 的 两 个 相 邻 单元 间 才 
可 能 产生 流 ; 若 单 元 i 不 是 汇集 单元 (中 心 单 


口 或 总 任务 数量 ， 商 为 平均 距离 。 对 于 人 口 
均等 分 区 , 目标 函数 值 表 示 区 域内 所 有 人 到 
达 其 中 心 点 的 平均 距离 。 约 束 条 件 (2) 保 证 每 
个 单元 必需 划分 在 一 个 分 区 中 。 约 束 条 件 (3) 
要 求 分 区 在 一 定 误差 内 均等 ， 式 中 0 = 
Diev pi/K 为 人 口 或 任务 量 平均 值 ，e 为 限 差 ， 
根据 实际 需求 取 值 ， 如 1%、5% 或 其 他 值 。 约 
束 条 件 (4) 要 求 每 个 分 区 数量 为 尺 。 约 束 (5)- 


元 )， 约 束 (7) 保 证 该 单元 至 少 产 生 一 个 单位 
的 流 。 约 束 (5)-(7) 将 使 每 一 个 分 区 中 的 产 流 
汇集 到 中 心 单 元 。 条 件 (9)-(1D 定 义 决 策 变 


三 | 


量 。 

通过 事先 确定 分 区 中 心 ( 流 汇 ) 单 元 可 将 
模型 简化 。 给 定 集 合 S={s1,s2...sSK}( Sk E 
Te Vk E 有) 表示 天 个 分 区 的 中 心 单 元 。 移 除 
决策 变量 y,， 以 上 模型 可 简化 如 下 : 


Min. Biey Dres Didis,Xis, / Liev Pi (12) 
S.T. YeesXn = LViEV (13) 
(1-£)Q < Yiev Dixis, < (1 + £)Q, Yk E€ S 14 
fije = Nn — K)xir Yİ EV,j EN ,kEsS (15) 
fiix S (n— K)xjp Vi EV,j EN, KES (16) 
Djeni fijk — bjen; fjik Z Xin, Vİ EV,kES,i +k (17) 
Xn = (0,1 Vi EV,kEsS (18) 
fij 2O,ViEV, JEN, RES (19) 


本 质 上 , 模型 (1)-(11) 相 当 于 增加 了 空间 
连续 约束 的 Location-Allocation 问题 模型 ; 
而 简化 模型 (12)-(19) 相 当 于 增加 了 空间 连续 
约束 的 Allocation 问题 。 与 原 问 题 相 比 ， 简 
化 问题 不 需要 确定 分 区 中 心 , 且 决 策 变量 数 
量 大 幅 减 少 。 也 应 带 注意 到 ， 容 量 约束 
Allocation 问题 仍 是 NP-Hard 问题 ， 简 化 模 
型 仍 具 有 很 高 的 计算 复杂 度 。 文 献 [37] 利 用 
K-medoids 算法 确定 分 区 中 心 ， 再 使 用 与 
(12)-(19) 接 近 的 模型 进行 分 区 ， 计 算 效 率 较 
高 。 

2 算法 设计 
2.1 基本 算法 


选择 迭代 局 部 搜索 (ILS) 算 法 作为 求解 
均等 分 区 问题 的 算法 框架 。ILS 算法 思路 简 
单 、 易 于 实现 ， 对 于 离散 优化 问题 行 之 有 效 
[38]。 该 算法 从 一 个 初始 解 开 始 ， 返 代 地 进 
行 扰动 和 局 部 搜索 。 局 部 搜索 容易 陷入 局 部 
最 优 , 对 当前 位 置 的 扰动 能 够 使 算法 脱离 局 
部 最 优 。ILS 算法 中 的 核心 部 件 包括 始 解 生 
成 、 局 部 搜索 和 当前 解 扰 动 ， 示 意 如 下 : 


Algorithm 1: Iterated local search 


1. so = GenerateInitialSolution(); 
2. s = LocalSearch(so); 

3. Repeat until ter. condition met: 
4. s' = Perturbation(s); 

5. s* = LocalSearch(s’); 


6. If f(s*) <f(s): s = 8*; 

7. Output s. 

有 多 个 方法 产生 初始 解 , 例如 随机 种 子 
生长 算法 和 加 权 K-medoids 算法 。 随 机 种 子 
生长 算法 原理 如 下 : 随机 选择 玉 个 空间 单元 
作为 分 区 中 心 , 按 距离 中 心 最 近 原 则 将 所 有 
单元 指派 给 分 区 中 心 ， 从 而 形成 初始 分 区 。 
加 权 K-medoids 算法 [37] 是 一 个 加 权 聚 类 算 
法 ， 随 机 选择 天 个 空间 单元 作为 分 区 中 心 ， 
迁 代 进行 加 权 指 派 和 中 心 点 更 新 ,直到 满足 
终止 条 件 。 这 两 种 算法 均 不 保证 分 区 连续 要 
求 和 分 区 均等 要 求 。 一 般 来 说 ， 两 种 算法 获 
得 分 区 结果 接近 空间 连续 , 仅仅 作 少 量 的 单 
元 调整 即 可 修复 分 区 为 连续 。 分 区 均等 条 件 
得 不 到 满足 时 ， 在 后 续 搜索 中 逐步 改进 ， 直 
到 满足 要 求 。 通常 的 做 法 是 在 目标 函数 中 加 
入 非 均 等 惩罚 项 ,， 随 惩罚 项 数值 降低 直至 为 
零 ， 分 区 均等 条 件 得 到 满足 。 

ILS 局 部 搜索 采用 分 区 边界 单元 移动 的 
方法 [39]， 包 括 每 次 移动 一 个 单元 的 一 单元 
移动 算 子 、 每 次 移动 两 个 单元 的 二 单元 移动 
算 子 ,以 及 每 次 移动 三 个 单元 的 三 单元 移动 
算 子 。 考 虑 到 每 次 移动 三 个 单元 的 计算 复杂 
度 很 高 ， 本 文 仅 采 用 前 两 个 算 子 。 

对 于 空间 均等 分 区 ， 扰 动 方法 很 多 。 例 
如 ， 破 环 若 干 分 区 、 破 环 知 干 连续 分 区 、 破 
环 一 个 连续 区 域 、 破 环 一 定 比 例 的 边界 单 
元 ， 然 后 进行 修复 。 也 可 以 直接 将 若干 边界 
单元 移动 到 相 邻 分 区 中 , 之 后 根据 需要 进行 
连续 性 修复 。 

2.2 算法 扩展 

为 提升 优化 性 能 , 对 ILS 算法 进行 三 方 
面 的 扩展 。 首 先 ， 改 造 单 解 ILS 算法 为 群 解 
ILS 算法 ， 其 次 ， 在 ILS 迭代 循环 中 ， 使 用 
变 邻 域 下 降 (VND) 搜 索取 代 简 单 的 局 部 搜 
索 ; 第 三 , 记录 ILS 搜索 过 程 中 产生 的 分 区 ， 
ILS 搜索 完成 后 使 用 SPP 模型 ， 进 一 步 选择 
更 优 的 方案 。 扩 展 后 的 算法 是 一 个 混合 算 
法 ， 流 程 如 下 : 

Algorithm 2: ILS-based hybrid algorithm 
参数 : 群 大 小 (psize)， 破 环 强度 
(strength)， 连 续 未 更 新 最 好 解 循环 数 
(mloops), SPP 模型 求解 时 间 (2). 


1. P = GenerateInitialSolutions(psize); 


2. pool=null; 

3. Sbes=Best(P); 

4. notImpr=0; 

5. While notlmpr < mloops: 

6. Select a solution s from P randomly; 
7. s'=Perturbation (s, strength); 

8. s"=VNDsearch(s’, /sops); 

9. s'=updateCenetrs(s”’) 

10. If fs*)<f(spes): Sbes™s*, notImpr=0; 
11. else: notImprt+=1, 

12. pool=UpdateAreaPool(pool, s*); 
13. P =UpdatePopulation (P, s*); 

14. s=SetPartitioning (pool, t); 

15. Output s. 


与 基于 单 解 的 搜索 算法 相 比 , 本 文 算法 
维护 一 组 解 。 首先 , 算法 步骤 (1) 生 成 一 组 初 
始 解 ; 其 次 ,在 每 一 次 迭代 ， 从 群 解 中 随机 
选择 一 个 解 进行 搜索 (步骤 6); 第 三 , 搜索 完 
成 后 ， 使 用 新 解 更 新 群 解 (步骤 13)。 群 解 更 
新 中 ， 优 先 考虑 解 的 目标 值 ， 其 次 考虑 解 的 
差异 程度 ， 保 持 解 之 间 有 一 定 的 差异 。 基 于 
群 解 的 ILS 算法 维护 一 组 具有 差异 度 的 精英 
解 ， 扩 大 了 解 空间 搜索 范围 ， 有 利于 改进 求 
解 质量 ;， 同 时， 算法 收敛 速度 通常 会 变 慢 ， 
计算 时 间 有 一 定 的 增加 。 

ILS 中 使 用 VND 搜索 能 够 使 每 轮 搜索 
达到 局 部 最 优 。VND 方法 反复 使 用 邻 域 搜 
索 算 子 进行 搜索 , 直到 所 有 算 子 均 不 能 提升 
当前 解 ， 从 而 获得 局 部 最 优 解 。ILS 中 频繁 
使 用 扰动 方法 ， 若 当前 解 不 是 局 部 最 优 ， 下 
轮 搜索 中 很 容易 被 破坏 , 减少 了 到 达 局 部 最 
优 的 机 会 。 因 此 ，VND 搜索 常用 于 ILS 算 
法 。 

算法 步骤 (13) 使 用 了 SPP 模型 。SPP 模 
型 求解 时 需要 先 构造 大 量 分 区 ， 再 通过 模 
型 选择 最 优 的 分 区 组 合 。ILS 搜索 过 程 中 ， 
会 产生 大 量 的 分 区 。 所 有 的 分 区 构成 集合 
Q。 每 一 个 分 区 i 具有 目标 值 O 和 单元 集 
A Ur SPP 模型 从 集合 8 中 选择 一 个 数量 
为 KK 的 子 集 ， 宪 盖 所 有 地 理 单元 集合 V, 
并 使 目标 值 最 小 。 定 义 2 为 8 的 子 集 ， 
0 = {ilie 0,j € Ui}， 是 一 个 包含 空间 单 
元 j 的 分 区 的 集合 ， 定 义 决策 变量 x 表明 
候选 服务 分 区 i 是 否 被 选中 。SPP 模型 如 
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Az 
Min. Vien OiXi (20) 
S.T. Diea; Xi 二 1, Vj EV (21) 
dicen Xi = K (22) 
x, = {0,1}, vie 0 (23) 

3 实验 验证 


3.1 实验 设计 


AN 


HR 


2 
S 


收集 5 个 区 域 的 地 理 数据 进行 均等 分 区 
实验 。 案 例 区 域 范 围 小 到 一 个 住宅 小 区 ， 大 
到 省 级 行政 区 ; 空间 单元 个 数 从 28 到 2145; 
既 有 城市 也 有 农村 地 区 。 案 例 区 形状 、 空 间 
单元 及 单元 属性 值 分 布 如 图 1 所 示 ,， 自 上 而 
下 从 左 到 右 区 域名 称 简 称 为 HD、ZY、GY、 
GY2 和 HN。 主 要 指标 统计 见 表 1。 


图 1 实验 区 示意 图 
Figure 1 The study areas 
表 1 实验 区 主要 特征 


Table 1 Main features of the study areas 


案例 区 域 HD GY ZY GY2 HN 
地 理 范围 居住 小 区 县 级 市 市 辖区 县 级 市 省 
空间 单元 楼 栋 行政 村 小 区 自然 村 乡镇 街道 
空间 单元 数量 28 297 324 1,276 2,145 
单元 属性 值 之 和 2,420 36,824 3,873 819,812 96,918,036 


对 五 个 区 域 进行 均等 化 分 。 考 虑 到 地 理 
空间 不 可 能 严格 均 分 ， 人 允许 分 区 有 一 定 的 误 
差 , 即 模型 中 的 参数 s。 实验 中 , 各 案例 区 e 取 
HA 5%。 因 案例 GY 和 GY2 中 空间 单元 属 
性 差异 很 大 , 分 区 数量 不 能 过 多 ， 否则 难以 
满足 均等 约束 。 

实验 计算 环境 为 HP 桌面 计算 机 ， 配 
置 Intel Core I7-6700 CPU 3.40-GHz 和 8GB 
WIF, Windows 10 操作 系统 , 安装 有 Python 
2.7、ArcGIS 10.4 和 IBM CPLEX 12.6 优化 
器 。 在 ArcGIS 中 准备 案例 数据 ， 包 括 空 间 


— 


单元 及 其 属性 pix 单元 邻接 关系 cy 和 距离 
dx 以 单元 中 心 点 为 基准 )。 本 文 算法 使 用 
Python 程序 设计 语言 编程 实现 。 为 提升 计算 
速度 ， 算 法 在 PyPy7 (https://www.pypy.org) 
环境 中 运行 。 
3.2 计算 结果 

针对 给 个 案例 区 域 和 分 区 数量 使 用 两 
种 方法 求解 。 第 一 ， 按 照 数 学 表达 式 (1)-(11) 
构建 案例 模型 ， 使 用 CPLEX12.6 优化 器 进 
行 求解 。 因 区 域 GY2 和 HN 规模 过 大 ， 仅 
完成 了 区 域 HD. GY 和 ZY 案例 的 模型 求 
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解 。 第 二 ， 使 用 2.2 节 算 法 求解 。 算 法 参数 
设置 如 下 : psize=10 ， strength=5% , 
mloops=100, SPP 模型 求解 时 间 限 制 二 100s。 
为 验证 算法 有 效 性 , 每 个 区 域 的 每 个 分 区 均 计 


案例 区 HD 分 区 结果 见 表 2, 分 区 示意 图 
见 图 2, KF, LB 和 UB 分 别 表 示 目 标 函 数 (1) 
的 下 界 和 上 界 ，MIPGap 为 CPLEX 求解 结 
A i le de 人 


算 10 次 ， 统 计 其 平均 目标 值 和 标准 差 ， 并 与 
模型 计算 提供 的 目标 值 下 界 作 比较 。 


因 该 区 域 规模 较 小 ， 
解 ， 混 合算 法 也 能 够 


表 2 区 域 HD 分 区 结果 


Table 2 Solution results from area HD 


模型 计算 较 快 获得 
快速 获得 最 优 解 。 


分 区 BR CPLEX 混合 算法 
量 2 LB UB MIPGap Time/s Gap Dev Time/s 
2 5% 0.099 0.099 Opt. 0.99 0.00% 0.00% 0.32 
3 5% 0.084 0.084 Opt. 1.20 0.00% 0.00% 0.44 
4 5% 0.071 0.071 i 1.14 0.00% 0.00% 0.52 
5 5% 0000 0.060 Opt 131 0.00% 0.00% _ 0.62 __ 5% 0.060 0.060 t. 1.31 0.00% 0.00% 0.62 


案例 区 GY 分 区 结果 见 表 3， 
意图 见 图 3。 可 以 看 出 ，CPLEX 优化 器 能 够 求 
解 空间 单元 数量 为 297 的 分 区 案例 ,文献 [5] 中 


图 2 区 域 HD 分 区 示意 图 
Fig. 2 Districting maps of area HD 
部 分 分 区 示 
其 接近 最 优 解 ; 
够 获得 近似 最 优 解 ; 


7 和 10 时 求解 质量 高 ， 能 够 获得 


求解 类 似 问 题 的 模型 ， 难 以 求解 36 空间 单元 。 IRE. 与 模型 精确 求解 相 比 ， 混 合算 法 效率 
的 分 区 案例 。 本 文 混合 算法 分 区 数量 为 2、3、 非常 高 ， 仅 需 6~10 秒 。 
表 3 区 域 GY 分 区 结果 
Table 3 Solution results from area GY 
分 区 BR CPLEX 混合 算法 
数量 E LB UB MIPGap Time/s Gap Dev Time/s 
2 5% 8.149 8.149 Opt. 6,098.57 0.00% 0.00% 6.20 
3 5% 6.085 6.085 Opt. 3,990.47 0.04% 0.06% 6.18 
4 5% 5.135 3.135 Opt. 11,145.66 1.19% 0.03% 7.51 
5 5% 4.616 4.616 Opt. 7,588.53 0.16% 0.07% 9.92 
6 5% 4.185 4.185 Opt. 11,303.26 0.24% 0.37% 6.67 
7 5% 3.751 3.751 Opt. 3,576.36 0.00% 0.00% 6.17 
8 5% 3.529 3.529 Opt. 5,390.92 0.22% 0.09% 6.92 
10 5% 3.120 3.120 Opt. 4,823.14 0.00% 0.01% 7.18 


最 优 解 或 极 
分 区 数量 为 5、6 和 8 时 ， 能 
分 区 数量 为 4 时 ,求解 质 
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图 3 区 域 GY 分 区 示意 图 
Fig. 3 Districting maps of area GY 
案例 区 ZY 分 区 结果 见 表 4， 部 分 分 区 结 。”” 仅 分 区 数量 为 20 时 ， 目 标 值 与 最 优 解 相 差 


果 示 意图 见 图 4。 可 以 看 出 ， 所 有 案例 模型 均 “0.29%。 另 外 ， 混 合算 法 稳定 性 较 好 ，10 次 计 
获得 最 优 解 ， 且 混合 算法 结果 有 逼近 了 最 优 解 ， 算 结果 的 标准 差 较 小 。 
表 4 区 域 ZY 分 区 结果 


Table 4 Solution results from area ZY 


分 区 BR CPLEX 混合 算法 
量 & LB UB MIPGap Times Gap Dev Time/s 


2 5% 0.979 0.979 Opt. 3,480.72 0.00% 0.00% 5.24 
3 5% 0.831 0.831 Opt. 3,796.54 0.00% 0.00% 3.19 
4 5% 0.717 0.717 Opt. 4,292.85 0.01% 0.04% 5.84 
5 5% 0.624 0.624 Opt. 4,758.05 0.05% 0.03% 5.80 
6 5% 0.562 0.562 Opt. 5,811.00 0.08% 0.02% 6.04 
T 5% 0.513 0.513 Opt. 5,927.70 0.01% 0.01% 6.00 
8 5% 0.477 = 0.477 Opt. 4,544.64 0.02% 0.04% 7.32 
10 5% 0.425 0.425 Opt. 4,973.77 0.02% 0.01% 7.03 
15 5% 0.351 0.351 Opt. 4,973.77 0.17% 0.10% 9.84 
20 5% 0.297 0.297 Opt. 5,397.47 0.29% 0.19% 15.54 


图 4 区 域 ZY 分 区 示意 图 
Fig. 4 Districting maps of area ZY 


案例 区 GY2 分 区 结果 见 表 5， 部 分 分 区 溢出。 混合 算法 求解 结果 稳定 性 较 好 ， 标 准 差 
示意 图 见 图 5。 案 例 区 GY2 和 HN 空间 单元 数 AEF 0.03%~0.11% 之 间 。 该 区 域 有 1276 个 空 
ERK, CPLEX 优化 器 求解 数学 模型 时 内 存 。” 间 单元 ， 能 够 在 2~5 分 钟 内 完成 分 区 计算 。 
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#5 区 域 GY2 分 区 结果 


Table 5 Solution results from area GY2 


分 区 数量 。 参数 e 最 小 平均 最 大 标准 差 。 Time/s 
5 5% 4.439 4.442 4.443 0.03% 115.78 
10 5% 3.108 3.109 3.110 0.03% 123.96 
15 5% 2.544 2.548 2.551 0.08% 268.82 
20 5% 2.174 2.177 2.181 0.11% 165.05 


2 had 


图 5 区 域 GY2 分 区 示意 图 
Fig. 5 Districting maps of area GY2 
案例 区 HN 分 区 结果 见 表 6， 部 分 分 区 示 。 ”标准 差 处 于 0.00%~0.18% 之 间 。 男 外 ,算法 需 
意图 见 图 6。 混 合算 法 求解 结果 稳定 性 较 好 ， 要 5~9 分 钟 完成 分 区 。 
表 6 区 域 HN 分 区 结果 


Table 6 Solution results from area HN 


分 区 数量 Be 最 小 平均 最 大 标准 差 。 Time/s 
5 5% 69.762 69.766 69.769 0.00% 302.46 
10 5% 46.548 46.560 46.586 0.02% 330.24 
15 5% 37.562 37.571 37.581 0.02% 546.56 
20 5% 31.781 31.795 31.810 0.03% 404.85 
25 5% 28.291 28.307 28.327 0.04% 539.99 
30 5% 25.854 25.864 25.879 0.03% 340.52 
35 5% 23.764 23.798 23.832 0.09% 323.62 
40 5% 22.244 22.296 22.333 0.13% 498.67 
45 5% 20.760 20.802 20.860 0.18% 388.39 
50 5% 19.541 19.590 19.612 0.12% 410.77 


图 6 区 域 HN 分 区 示意 图 
Fig. 6 Districting maps of area HN 
3.3 应 急 服 务 应 用 分 外 ， 均 等 分 区 还 能 够 广泛 应 用 于 应 急 管 理 、 
除 应 用 于 经 典 的 选区 、 销 售 区 和 巡逻 区 划 ”城市 管理 责任 区 划分 、 区 域 公 共 服 务 等 。 假定 


在 新 冠 病毒 流行 期 间 , 疾 控 部 门 启 动 某 区 域 全 
员 核 酸 检 测 , 派出 若干 小 组 完成 测试 。 均 等 分 
区 结果 可 满足 检测 基本 要 求 : 每 个 组 任务 均 
等 ， 能 快速 完成 检测 任务 ; 每 个 服务 站 服务 的 
居民 尽 可 能 集中 , 从 而 使 居民 到 达 服 务 站 步行 
距离 最 短 ; 避免 聚集 ， 即 居民 仅 和 近邻 居民 一 
起 排队 ， 避 人 免 与 居住 较 远 的 居民 聚集 在 一 起 ; 
方便 管理 ， 每 个 服务 点 的 服务 区 域 空间 连续 ， 
便于 发 布 通知 、 上 门 服务 等 。 案例 HD 适用 于 
居民 小 区 全 员 核 酸 检测 .在 图 2 中 各 分 区 中 心 
位 置 或 中 心 附近 设置 服务 点 , 即 可 较 好 地 完成 
服务 点 与 服务 区 布局 规划 。 

均等 分 区 也 适用 于 应 急 物资 配送 ,假定 在 
新 冠 病 毒 流行 期 间 ， 某 区 域 实施 封闭 管理 。 为 
保障 医药 或 生活 物资 供应 , 应 急 指挥 部 门 根据 
本 区 域 参 与 物资 供应 商家 及 自愿 者 情况 , 将 其 
划分 为 若干 物资 配送 小 组 , 每 个 小 组 负责 一 个 
区 域 的 物资 供应 ,均等 分 区 可 以 保证 每 个 小 组 
任务 量 接 近 ， 分 区 形状 紧 赎 能 够 减少 交通 时 
间 。 例 如 ， 区 域 ZY 和 GY 分 区 可 作为 城市 和 
农村 地 区 物资 供应 服务 区 划分 。 因 计算 速度 
TR, 指挥 部 门 可 以 根据 参与 商家 和 自愿 者 的 动 
态 变 化 不 断 调整 服务 区 ， 优 化 物资 配送 服务 。 
4 讨论 与 结论 

针对 均等 分 区 问题 , 本文 构 造 了 一 个 混 
合 整 型 线性 规划 模型 ， 并 设计 了 一 个 基于 从 
代 局 部 搜索 算法 的 混合 算法 。 选择 5 个 区 域 
对 算法 进行 测试 ， 结果 发 现 : 中 本 文 模型 能 
够 求解 空间 单元 数量 达到 324 的 案例 ; QA 
文 设计 的 混合 算法 优化 性 能 优异 ， 鲁 棒 性 
强 ， 计 算 效 率 较 高 。 

本 文 在 三 个 方面 区 别 于 现 有 模型 : OÑ 
过 给 定 分 区 中 心 单元 , 使 分 区 连续 约束 条 件 
得 到 简化 。 与 现 有 模型 相 比 ， 该 模型 不 需要 
引入 汇流 单元 决策 变量 [4][$]。 文 献 [外 ] 均 质 
分 区 问题 模型 仅 能 够 求解 空间 单元 数 为 
16~49 的 部 分 案例 ,文献 [5] 政 治 分 区 问题 仅 
能 够 求解 空间 单元 数 为 25~36 的 部 分 案例 ， 
而 本 文 模型 可 用 于 求解 单元 数量 达 324 的 较 
大 规模 案例 。@ 通 过 给 定 分 区 中 心 单元 ， 易 
于 建立 紧凑 度 目标 函数 。 分 区 形状 紧凑 度 的 
定义 方式 很 多 , 本 文 使 用 加 权 旅 行距 离 作为 
紧凑 度 目 标 ， 物 理 意 义 明确 ， 易 于 构造 线性 
函数 , 也 避免 了 文献 [4] 引 入 额外 的 虚拟 决策 


变量 。 包 本 文 模型 本 质 上 是 容量 约束 
Location-Allocation 模型 和 空间 连续 约束 的 
组 合 。 模 型 求解 结果 提供 中 心 区 位 ， 从 而 使 
均等 分 区 问题 能 够 应 用 于 多 种 服务 规划 , 特 
别 是 对 于 服务 设施 位 置 条 件 要 求 不 高 的 问 
题 。 例 如 ， 新 冠 疫情 应 急 管 理 中 ， 可 以 基于 
均等 分 区 进行 医药 配送 、 生 活 物 资 配送 和 核 
酸 检查 等 服务 规划 ， 较 好 地 满足 规划 要 求 。 

本 文 设计 的 混合 算法 是 对 ILS 算 法 的 扩 
E: 基于 群 解 搜索 、VND 搜索 及 SPP 模型 
改进 。 虽 然 案例 测试 表明 算法 性 能 优异 ， 但 
仍 需 进一步 研究 。 一 是 构造 更 多 的 案例 对 算 
法 进行 充分 测试 ， 二 是 进行 算法 分 析 ， 引 入 
学 习 机 制 ， 自 适应 地 调整 算法 参数 。 
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The Equal Districting Problem: Model, Algorithm and Applications 
KONG Yunfeng! 

1 Key Laboratory of Geospatial Technology for the Middle and Lower Yellow River Regions, Ministry of Education, Henan University, 
Abstract: Districting problems have been widely applied in geography, economics, environmental 
science, politics, business, public service and many other areas. The equal districting problem (EDP) 
arises in applications such as political redistricting, police patrol area delineation, sales territory design 
and some service area design. The important criteria for these problems are district equality, contiguity 
and compactness. A mixed integer linear programming (MILP) model and a hybrid algorithm are 
proposed for the EDP. The hybrid algorithm is designed by extending iterative local search (ILS) 
algorithm with three schemes: population-based ILS, variable neighborhood descent (VND) local search, 
and set partitioning. The performance of the algorithm was tested on five areas. Experimentation showed 
that the instances could be solved effectively and efficiently. The potential applications of the EDP in 
emergency services are also discussed. 

Key words: equal districting problem; mixed integer linear programming; hybrid algorithm; emergency 


service 
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